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Abstract 

In this paper, we present the controllability properties of Keplerian motion controlled 
by low-thrust control systems. The low-thrust control system, compared with high or even 
impulsive control system, provide a fuel-efficient means to control the Keplerian motion of a 
satellite in restricted two-body problem. We obtain that, for any positive value of maximum 
thrust, the motion is controllable for orbital transfer problems. For two other typical problems: 
de-orbit problem and orbital insertion problem, which have state constraints, the motion is 
controllable if and only if the maximum thrust is bigger than a limiting value. Finally, two 
numerical examples are given to show the numerical method to compute the limiting value. 
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1 Introduction 


In classical mechanics, the determination of the motion of two celestial bodies, which interact only 
with each other, is the typical two-body problem. If one body is light enough, the uncontrolled 
motion of the light body around a heavy body is a restricted two-body problem, and the motion 
is well-known as Keplerian motion. A common example is the artificial bodies, i.e., spacecrafts or 
satellites, moving around the Earth. Once the atmospheric effects are negligible and the Earth 
overwhelmingly dominates the gravitational influence, a satellite moves stably on a periodic orbit 
if the mechanical energy of the satellite is negative. As an increasing number of artificial satellites 
have been launched into space around the Earth or even into deeper space since the mid of last 
century, an important problem arises in astronautics, that is to control the Keplerian motion of a 
satellite to transfer between different orbits to achieve desired mission requirements. 

The control of an artificial satellite is generally performed by system propulsion, expelling mass 
in a high speed to generate an opposite reaction force according to Newton’s third law of motion. 
Up to now, there are already several types of propulsion systems available, including chemical 
propulsion systems and electric propulsion systems. Though chemical propulsion systems are able 
to provide much higher thrust, electric propulsion systems have the the potential for a much higher 
specific impulse than is available from chemical ones, resulting in a lower fuel consumption and 
thus a longer satellite lifetime for a given propellant mass. On the one hand, the electric propulsion 
systems provide a fuel-efficient means to control the motion of a satellite; On the other hand, the 
fact that the possible maximum thrust, which the electric propulsion systems can provide, is very 
low results that the transfer time is exponentially long. Hence, the optimization of transfer time 
has be studied in Refs. [4, 5, 7]. In addition, the strong-local optimality of such problems was 
studied in Ref. [16]. 

Though the low-thrust control systems provide an fuel-efficient means to control the motion of 
a satellite, the problem that whether or not it has the ability to move a satellite from one point 
to another one arises. This is actually a controllability property, which is a prerequisite to analyze 
mission feasibility during designing a space mission or designing an optimal trajectory. Restricting 
the mechanical energy of a satellite into negative region without any other state constraints, the 
controlled motion is called orbital transfer problem (OTP), and the controllability for OTP was 
derived in Ref. [5, 7] to show that there exists admissible controlled trajectories for every OTP if the 
maximum thrust is positive. In the current paper, the controllability for OTP is established using 
alternative techniques from geometric control, (cf. Refs. [9, 12]). Taking into account the state 
constraint that the radius of a satellite is larger than the radius of the surface of the atmosphere 
around the Earth, the orbital insertion problem (OIP) and de-orbit problem (DOP) are defined 
in this paper. Some controllability properties for OIP and DOP are then addressed and we show 
that there exist admissible controlled trajectories for OIPs and DOPs if and only if the maximum 
thrust is bigger than a specific value (depending on the initial point or final point). 

The organization of the paper is the following. In Section 2, we recall the basic properties of the 
dynamics of the motion of a satellite around the Earth, and basic notations and definitions are given 
which are crucial for analysis of controllability properties. In Section 3, the controllability for OTPs 
is represented firstly by using geometric control technology in Ref. [9]. Then, the controllability 
properties for OIPs and DOPs are derived in Subsections 3.2 and 3.3, respectively. In Section 4, 
two numerical examples are given to show the development in this paper. Einally, a conclusion is 
given in Section 5. 

2 Notations and definitions 

2.1 Dynamics 

Consider a satellite as a mass point moving around the Earth, its state in a geocentric inertial 
cartesian coordinate (GICC), illustrated by Eigure 1, consists of its position vector r G 
velocity vector t; G and mass m G Then, the dynamics for the movement of the satellite 


2 


for positive times can be written as: 




sat 


'r{t) = v{t), 

: < v{t) = - ||r(t)||3»'(^) + 

m{t) = -P II T{t) II, 


( 1 ) 


where /i > 0 is the gravitational constant, /3 > 0 is a scalar constant determined by the specific 
impulse of the low-thrust control system equipped on the satellite, || • || denotes the Euclidean 
norm and the thrust (or control) vector r G takes values in the admissible set 

T(Tniax) = {t e I || T ||< Tmax}, (2) 

where Tmax is a positive constant. 



Figure 1: Geocentric inertial cartesian coordinate. 

li X = M^\{0} X and x = (r, v), we define two vector fields /q, fi on A' by 

fo{x)=( _l), (3) 

f,{x)=(^^^y ( 4 ) 

with denotes the set of 6 x 3 matrices wih real entries and Is is the identity matrix of M^. 

Let Bs be the closed ball in centered at the origin and of radius 5 > 0. For every 5 > 0, we 
consider the control-afhne system given by 

: x{t) = fo{x{t)) + f^{x{t))u{t), (5) 

where the control vector n G takes values in We will use in this paper the vector field point 

of view of Refs. [11, 13, 14]. For every point x G X and every u G we denote by 

f:XxB^^T^X, {x,u) ^ f{x,u) = f^ix) X fi{x)u, (6) 

where /q and fi are referred to as the drift vector field and the control vector fields respectively. 
Note that trajectories of F^ starting at any Xq G X and measurable u : M+ ^ B^ are well-defined 
on an open interval of M+ containing 0, which depends in general on xq and n(-). 

2.2 Study of the drift vector field in A" 

In this paragraph, we recall the main properties of the drift vector field /q. For every x G X, we 
use 7 a 5 to denote the restriction to M+ of the maximal trajectory of /q starting at cc, i.e. is 
defined on some interval [0,t/(cc)) where tf{x) < oc. Then the follwing holds true. 

Property 1 (First integrals [2, 6]). For every x e X, if ^x{t) = (f(t),'r(t)) on [0,t/(x)), the 
quantities 

h{t) = f{t) X v{t), (7) 

vih) 

L{t) = v{t) X h - (8) 
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(9) 


Ht) II ’ 


are constant along 7 ^ and the corresponding constant values are the angular momentum vector 
h the Laplace vector L G and the mechanical energy of a unit mass E which is the 

sum of the relative kinetic energy || v{t) |p /2 and the potential energy —fi/ || r{t) ||. 

As a consequence of Eq.(7) and Eq.( 8 ), we have the following two properties. 

Property 2 (Straight line [ 2 , 6 ]). Let x ^ X with h = 0, i.e. r and v are colinear. Then the 
trajectory 'Is a straight line and tf{x) is either finite or infinite. 

Property 3 (Conic section [ 2 , 6 ]). Let x e X with h ^ 0 i.e. r and v are not colinear. Then, 
tf{x) = and the trajectory 'is a periodic trajectory with locus defining a conic section lying 
in a two-dimensional plane perpendicular to h called the orbital plane. 


Let 


X = {(r, v) G A I r X v = / 0}. (10) 

Define on X the function e : x ^\\ L || /pi. Along every trajectory of /q starting at cc G A, one 
gets, after multiplying Eq.( 8 ) by f(t), that 


II 


II h 11^ 

/i(l + e{x) cos 0{t, x)) ’ 


( 11 ) 


where the angle 0(t,x) is defined by cos 0{t,x) = Note that the previous formula holds 

true if L = 0 since in that case e{x) cos 0{t, x) is equal to zero and the orbit is a circle. 



Eigure 2 : The orientation of a 2 -dimensional orbital plane in GICC and the geometric shape and 
orientation of an elliptic orbit on the orbital plane. 

Notice from Eq.(ll) that an orbit {r{t),v{t)) = 7x(t) with x G X on M+ is a parabola if 
e{x) = 1 and a hyperbola if e{x) > 1. Put a satellite on a parabolic or a hyperbolic orbit without 
no control, then it can escape to infinity lim r{t) = + 00 . Thus, parabolic and hyperbolic orbits 

t^+oo 

are generally used for a satellite to escape from the gravitational attraction of Earth. Eor every 
point X G A, ifO < e{x) < 1, the orbit 7x(^) on is an ellipse, whose orientation is illustrated 
in Eigure 2 . Moreover, it is easy to deduce the following characterization of elliptic orbits. 

Property 4. Given every point x G X, the mechanical energy E is negative if and only ife{x) < 1. 
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Thus, let us define the set 


V = {xeX\E<{)}, (12) 

then for every point {r^v) G 7^, the associated orbit on IR+ is periodic and the set V is called 
the periodic region in X. 

Definition 1 (Smallest period tp). Given every point x eV, we denote by 

tp I E —y M, X I —y tp{x^^ 
the smallest period of the orbit on M+. 

According to Eq.(ll), for every point cc G 7^, if e{x) 7^ 0, the associated orbit jx on [0?^p(^)] has 
its perigee point and apogee point at 6>(t, x) = 0 and tt, respectively. Thus, let 

Tp : P ^ E, T„ix) = TriW' 

r„:P^K,rp(x) = (M) 

we say rp{x) and ra{x) are the perigee and apogee distances of the orbit jx on [0, tp{x)] if e{x) 7^ 0. 
Note that ra{x) = rp{x) if and only if e{x) = 0, which corresponds to a circular orbit. 

Property 5 (Minimum radius and maximum radius). Given every periodie orbit (f(t),'h(t)) = 
7aj(t) on [0,tp(cc)] in V, we have rp{x) <|| r{t) ||< ra{x) on [0,tp(cc)]. Thus, the perigee dis- 
tanee rp(x) and apogee distanee ra{x) are the minimum radius and maximum radius of the orbit 
{f(t),v{t)) on [0,tp{x)]. 

2.3 Admissible controlled trajectory of Egat 

For every initial point = (cc^, mi) ^ X xM^ and measurable control function t(-) taking values in 
T(Tmax) with Tmax > 0, let t/ G bc the maximum time such that the corresponding trajectory 
r(t,T,yJ of Fsat lies in A X M^, i.e., if r(t,T,yJ = (x(t), m(t)), then x{t) G X and m{t) > 0 on 
[0,tf)' We use Xr to denote [0, tj). We say r(t, r, m^) on Tp is the controlled trajectory of Fgat 
starting from yi and associated with t(-). 

Remark 1. For every point yi = {xi,mi) G 7^ x let {x(t),m(t)) = r(t,0,2/J onXr, we have 
that x(t) = 7xi(t) and m(t) = mi for every t >0, i.e., tf = oo. 

Definition 2 (Controlled Keplerian motion). Given every initial point yi = {xi,mi) G 7^ x R'f 
and measurable eontrol funetion t(-) taking values in T{r^ax) with Tmax > 0; the eorresponding 
trajeetory r(t, t(-), m^) o/Fgat is ealled a eontrolled Keplerian motion. 

Let Tc > 0 and Mq > 0 denote the radius of the surface of atmosphere around the Earth and 
the mass of a satellite without any fuel, respectively, then given every point {x,m) G V x on 
the trajectories of Keplerian motions, it is required that || r ||> Tc and m > Mq. 

Definition 3 (Admissible region). We define the set 

A= {x = {r,v) gV \ ||r||>rc}, (15) 

the admissible region in V for Keplerian motion and/or eontrolled Keplerian motion. 

Definition 4 (Admissible controlled trajectory). Given every Mq > 0, we say the eontrolled 
trajeetory (x{t),m(t)) = V(t,T,Xi,mi) 0/ Fgat on some finite intervals [0,tf] C Xp with initial 
eondition {xi^m/ G Ax is an admissible eontrolled trajeetory if {x{t)) G A and m{t) > Mq 
fort G [0,t/]. 

For every time interval [0,t/] C Xp, since m{t) < 0, it follows m{tf) < m{t) on [0,t/]. Thus, the 
inequality m{t) > Mq can be ensured by m{tf) > Mq. 
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2.4 Controlled problems in A 

For X e A, let (f(t),'r(t)) = jx(^) on we have that the inequality || r(t) ||> Vc is satisfied on 
M+ if rp(x) > Tc. Thus, we define the set: 

p+ = {(r, v) eV : rp{x) > Vc}. (16) 

It is immediate to see that the periodic uncontrolled trajectory 7(t, x) starting at any x G 
remains in . 

Let 


p = {x = {r,v)eV\ \\r\\>rc,rp{x)<rc<ra{x)}. (17) 

Then, for every point x G V~, there exists an interval [^1,^2] G [0,tp(cc)] such that || r{t) ||< Tc for 
t G [^1,^2]. Thus, placing a satellite on a point x G V~, it can move out of the admissible region 

Definition 5 (Stable periodic region and unstable periodic region V~ in A). We say that the 
two sets and V~ are the stable and unstable periodie regions, respeetively. 

All the satellites periodically moving around the Earth are located in the stable periodic region 
. In order to fulfill observation or other mission requirements, a satellite is controlled to move 
from one point Xi in to another point Xf in by its control system. 

Definition 6 {Orbital Transfer Problem (OTP)). We say that the problem of eontrolling a satellite 
from a point Xi in to another point Xf in is the orbital transfer problem, see the first figure 
of Fig. 3a. 

For a typical space mission, in order to place a satellite into a stable orbit in a rocket is used 
to carry the satellite from the surface of the Earth to a point Xi in V ~, at which the rocket and the 
satellite are separated. From this moment on, the satellite is controlled by its own control system 
to be inserted into a stable orbit in . 

Definition 7 {Orbital Insertion Problem (OIP)). We say that the problem of eontrolling a satellite 
from an initial point Xi G V~ to a final point Xf G is the orbit insertion problem, see the third 
figure of Fig. 3b. 

After a satellite in the stable region finishes its mission, it should be decelerated to return to 
the unstable region V~. Then, the satellite will coast into atmosphere such that the aerodynamic 
pressure will act as a control to control the satellite to fly to landing sites. 

Definition 8 {De-Orbit Problem (DOP)). We say that the problem of eontrolling a satellite from 
an initial point Xi G to a final point Xf ^ V~ is the de-orbit problem, see the seeond figure of 
Fig. 3e. 


3 Controllability 

According to the definition for controlled Keplerian motion in Definition 2, the controllability of 
Keplerian motion deals with the existence of admissible controlled trajectories for OTP, OIP, and 
DOP. 

Definition 9 (Controllability for OTP). We say that the system Egat is eontrollable for OTP if 
there exists Tmax > 0 50 that, for every initial mass > 0 and every initial and final points 
{xi,Xf) G {V^)‘^, there exists a time tf G Tr cind an admissible eontrolled trajeetory {x{t),m{t)) = 
T{t,T,Xi,mi) o/Egat on [0,t/] inA^Mfi^sueh that x{tf) = Xf. 

Definition 10 (Controllability for OIP and DOP). We say that the system Egat is eontrollable for 
OIP (DOP respeetively) from any point Xi G V~ (xi G respeetively) if for every initial mass 
mi > 0 there exists Tmax > 0 50 that, for every final point Xf G (^f ^ respeetively), there 

exists a time tf G Tr ci'^d an admissible eontrolled trajeetory {x{t),m{t)) = V{t,T,Xi,mi) 0/Egat 
on [0,tf] in Ax sueh that x{tf) = Xf. 
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(a) OTP 


(b) OIP 

Figure 3: OTP, OIP, and DOP. 


(c) DOP 


For every initial point Xi ^ X and 5 > 0, we use V^(t^u{t)^Xi) to denote the trajectory of 
in Eq.(5) associated with a measurable control u{-) : [0,f/] ^ Be and we define tf G as 
the maximum time such that Te{t^u{t)^Xi) lies in T' on [0,t/). Set X = [0,f/). We refer to 
Te{t, u{t), Xi) as the controlled trajectory of starting from Xi and corresponding to the control 

m(-). 

Remark 2. Since ^x{t) = Te{t^0,x) onX, the uncontrolled trajectory Te{t,0,x) is periodic on 
if X G V. 

Lemma 1. Fix 5 > 0 and = {xi^irii) e X x R'f. Then, given every measurable control 
u(’) : [0,t/] ^ Bey if Tmax > ^TUi, then there exists Mq > 0 and an admissible controlled trajectory 
{x{t),m{t)) = r{t,T,Xi,mi) o/Ssat on [0,t/] in A x [Mo,mi] such that re{t,u,Xi) = x{t) for 
every t G [0,tj] and rri{tf) > Mq. 

Proof Since m{t) < mi for each time t G [0,t/] and Tmax ^ it follows that the thrust vector 
t(-) on [0,t/] can take values in the set T in Eq.(2) such that T{t)/m{t) = u{t) forevery time t G 
[0, t/]. Thus, let {x{t),m{t)) = r(t, r, Xi,mi), we have Te{t, u, Xi) = x{t) for every time t G [0, t/]. 
Since along the trajectory {x{t),m{t)) = T{t,T,Xi,mi) on [0,t/], we have u{t) = T{t)/m{t), which 
implies that m{t) = —/3 || u{t) || m{t). Thus, we obtain 

mit) = 

> ( 18 ) 

Let Mo := > 0. Then m{tf) > Mq and the lemma is proved. □ 

In order to study controllability, it is necessary to first show that the admissible region ^ is a 
connected subset of V. 

Lemma 2 (Connectedness of A). The admissible region A is an arc-connected subset of V, i.e., 
for every initial point Xi ^ A and every final point Xf ^ A, there exists a continuous path f : 
[0, 1] ^ A, A i-G x{X) such that x{0) = Xi, and x{l) = Xf. 

Proof We use the MEOE coordinates (cf. Definition 15) to prove the result, i.e., it is enough to 
show that Z is arc-connected. Let us choose two point 2^^ and 2;/ in ^ given by 

hx^y hy^,li) , Zj {^Pf , Cx ^ y Cy ^ , hx f y hy ^ ,l , 

with Xi = x{zi) and Xf = x{zf). We thus define the path z : [0,1] ^ ^ by 

;2(A) = (P(A),e,(A),e,(A),h,(A),h,(A),/(A)) 

where 

P(A) = [(1 - X)ri + Ar/][1 + ex{X) cos(/(A)) + e^(A) sin(/(A))], 

^ai(A) — (1 A)e3;;^ T Xcxj^y CyiX^ — (1 A)ey^ T Xcy^, 

hx{X) = (1 — X)hx^ + Xhxfy hy{X) = (1 — X)hh. + (1 — A)ea,., 
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where Vf =|| Vf || and =|| ||. Note that ex{X)^ + < 1 for each A G [0,1]. Let 

g{X) = (P(A), ea,(A), ey(A), ha,(A), hy(A),/(A)) on [0,1], we then have that ^(0) = 2:^ and ^(1) = 
Zf. Consider the continuous function x{X) = {r{z{X)),v{z{X))) for A G [0,1]. It follows that 
cc(0) = Xi and x{l) = Xf. It is immediate that x{X) G ^ for A G [0,1]. Finally, since P(A) = 
(r(A) X v{X)Y/fi > 0 and 0 < e = y/e^^Ap^Fe^^A)^ < 1, we have r(A) x v{X) F 0 and E{X) = 
- ||t»(F)|| < 0 on [0,1]. This proves the lemma. □ 

We also need the following lemma. 

Lemma 3 (Connectedness of P+). The set P+ is a eonneeted subset of A. 

Proof Given every two points Xi G and cc/ G P+, using the same technique as in the proof of 
Lemma 2^ let x{X) = {r{z{X))^v{z{X))) on [0,1], but we rewrite P(A) in the following form. 


P(A) = ((1 - A)rp, + Arp,)(l + fTPf+LAf), 

where r^. = rp{xi) and = rp{xf). Then, we have 

P(A) 


rp{x{X)) = 


1 + y ex{\)^ + ey{X)^ 

Thus, x{') takes values in P+ and this proves the lemma. 


= (1 - X)rp. + Xrp > Tc. 


□ 


3.1 Controllability for OTP 

In this subsection, we first give a controllability property of for OTP, then, according to Lemma 
1, we will establish the controllability of Sgat for OTP. 

Definition 11. For every eontrolled trajeetory x{') = T^^'^u^Xi) of (where 5 > 0, 'u(') : 
[0, tj] ^ Bs is measurable and Xi G A), we define 

S*(x) : X{t) = A{t)X{t) + B{t)u{t), 

the linearised system along x(-) ofE^ on [0,^/], where 

Mt) = B{t) = f^{x{t),u{t)), 


on [0,t/]. 

We first have a result of local controllability for the systems S^’s around the periodic trajectories 
of the drift vector field. 


Lemma 4. Let x G P+. Then, for every p > 0, there exists a > 0 sueh that the following properties 
hold: Ba{x) C P+ and, for every x G Ba{x), there exists a eontrolled trajeetory V^(t,u{t),x) of 
Eg sueh that 

re(0, u, x) = X, r^{tp{x),u, x) = X, 


and 

fort e [0,tp(S)]. 


II r^{t,u,x) -r^{t,o,x) ||< p, 


Proof According to Theorem 7 of Chapter 3 in Ref. [12], it suffices to prove the controllability of the 
linearized system E*{T^{t,0,x)) along the periodic trajectory T^{t,0,x) on the interval [0,tp{x)]. 
Then, the latter controllability would follow, according to Corollary 3.5.18 of Chapter 3 in Ref. [12], 
by the following rank condition: there exists a time r G [0, tp{x)] and a nonnegative integer k such 
that the rank of the matrix [Bo{r), Pi(r), • • • , Bk{r)] equals 6, where = A{t)Bi{t) — ■^Bi{t) 

for i = 1,2, • • •, and Po(F = therefore amounts to compute some The explicit 

expressions for matrices A and B in terms of x are 


A = f., 


II ngiS -h O II ||5 f f 


h 

0 


and B = 


0 

h 
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Since Bo{t) = B(t)^ it follows that 


B^{t) = A{t)Bo{t) - j^Boit) 

= A{t)Bo{t)=[h 0]^. 


Thus, we have that the rank of the matrix [Bo{t)^ Bi{t)] = 
t G [0,tp(x)], proving the lemma. 


is equal to 6 for every time 

□ 


Proposition 1. For 5 > 0, the control system is controllable for OTP within , i.e., for 
every initial point Xi G and final point Xf G , there exists a controlled trajectory re(t, n, Xi) 
ofTi^ in on a finite interval [0,t/] CX such that Ts{tf^u^Xi) = Xf. 


Proof Since the subset is path-connected as is shown by Proposition 3, it follows that any 
two different points Xi and Xf in are connected by a path x : [0,1] ^ P+, A i-^ x{X) such 
that cc(0) = Xi and x{l) = Xf. By compactness of the support of cc(-) in P+ there exists, for 
every cr > 0, a finite sequence of points Xq^Xi, • • • ,Xn^ on the support of x{X) so that Xq = 

Xn = Xf and G Bcj{xj)^ for j = 0,1, • • • — 1. According to Lemma 4, for every p > 0, 

there exists a > 0 small enough and a finite sequence of points Xq^Xi^ - • • , xat as above such that, 
for j = 0,1, • • • , — 1, Xjj^i G B(j{xj) and one has a controlled trajectory T^it^Uj^Xj) on the 

interval [0^tp{xj)] such that 


Tg (0, Uj , Xj ) — Xj , Tg (fp {xj ), Uj , Xj ) — x , 


and 

II V^{t,Uj{t),Xj) - r^{t,0,Xj) ||< p, for t G [0,tp{xj)]. 

For cr > 0, let Wj C P+ be an open neighborhood of Xj such that B(j{xj) C Wj for j = 0,1, • • • 
and set V^. = {r£(t,0,Wj), t > 0}. For p > 0 small enough, the open set is included in 
P+. By concatenating the T^{',Uj,Xj) for j = 0,1, • • • , A" — 1, the initial point Xi can be steered 
to Xf, proving the proposition. □ 

According to Lemma 1, and recalling the definition of controllability for OTPs in Definition 6, we 
obtain the following result of controllability: 

Corollary 1. For every /a > 0, /3 > 0, Tmax > 0, the system Ugat is controllable for OTP. 

Note that P+ C A, so the system Ugat is controllable for OTPs within A no matter what value of 
^max the low-thrust control system can provide if the satellite takes high enough percent of total 
fuel, i.e., {mi — mf)lmi > 0 is big enough, which makes senses in engineering for electric thrust 
systems whose maximum thrust r^iax is very small. 

3.2 Controllability for OIP 

We provide next a controllability criterium for OIP. 

Lemma 5. Assume that, for every point {xi,mi) G V~ x there exists r > 0 and a positive 
time f G Xr ci'^d a control t(-) : [0,t] ^ such that along the controlled trajectory {x,m{t)) = 

T{t,T{t),Xi,mi) on [0,t], we have x{t) ^ A on [0,t], m{t) > 0, and rp{x{t)) > Vc- Then, the 
system Ugat is controllable for OIP from {xi,mi). 

Proof Note that the assumption implies that there exists a control t(') : [0,t] ^ such that 
the admissible controlled trajectory {x{t),m{t)) = T{t,T,Xi,mi) in A x [m{t,mi] on [0,t] steers 
{xi,mi) in V~ x to some {x{t),m{t)) in P+ x After arriving at x{t) in P+, according to 
Proposition 1, it follows that there exists an Mq G (0,m(t)], a finite time G Xr, and a control 
t(-) : [0,t/] ^ T(r) such that along the controlled trajectory {x{t),m{t)) = r(t, r, x(r), m(r)) on 
[0,t/], we have x{t) G Al on [0,tj], x{tf) = Xf, and m{tf) > Mq. □ 
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One cannot have controllability for OIP for every value of Tmax > 0- Indeed, pick a point (cc^, m^) in 
V~ X For every control t(-) taking values in T(Tmax )5 the corresponding controlled trajectory 
(cc('), m(')) = r(', T, m^) converges to r(', 0, m^) on [0,tp(cci)] as Tmax tends to zero. Then, 

since rp{xi) < rc^ there exists t G [0,tp(cc^)] such that || r{t) ||< Tc implying that (cc('),m(')) = 
r(', T, cci, m^) is not admissible for every control t(') taking values in T(rniax)- For large values 
of Tmax: wc cau stccr (xi^nfii) G V~ X to (xf^rrif) G 7^^ x as described in the following 
lemma. 

Lemma 6. For every /3,/i > 0 and point = {xi^rui) in V~ x ^/^ere exists Tmax sueh that 
the following holds: 

1) if r > Tmax; there exists a eontrol t(-) taking values in T{r) and a positive time t G Xp sueh 

that, along the eontrolled trajeetory {x{t),m{t)) = r(t, r, 7/^, m^) on [0,t], we have x{t) G A 
on [0,f\, m{f) > 0; and rp{x{f)) > Vc', 

2) if r < Tmax; for every eontrol t(-) taking values in T{r), the eontrolled trajeetory (x(-),m(-)) = 

,mi) does not reaeh x 


Proof. At the light of the remark preceding the lemma, it is enough to find a value of t and a 
control t(-) taking values in T(t) steering = {xi,mi) to some point in x along an 
admissible trajectory. We proceed as follows. Let Cq = which belongs to (0,\/2p^). 

Choose now the function v{') defined on some time interval [0,T] (with T to be fixed later) as 
follows: i;(0) = Vi, r{t) = v{t) and 


v{t) 


(2r(i))V2 ^ *||r(t)|| ^ ^\\r{t)\\J 


te[0,T]. 


Here r{t)-^ denotes a continuous choice of vector perpendicular to r{t) in the 2D plane spanned by 
Vi and Vi. (Implicitely, we assume with no loss of generality that Vi 0 with > 0 and bi ^ 0.) 
For simplicity, we assume next that Oi = bi = 1. The curve (r(•),!;(•)) defined previously can be 
explicitely integrated using polar coordinates for r(-) = r(-) exp(i^(-)). One gets that, on [0,T], 

/ \ 1/2 


After integration, one has for t G [0,T], 

r{t) = (rf " + 0{t) = 0iF2 ln(rf" + 3Cot/2)/3. 

One also checks that h{t) = with a constant vector of unit norm parallel to x Vi 

and L{t) = {CqI 2 — y) ||^|^]|| • One deduces that rp{t) = for some positive constant Ci. 

One thus fixes T so that rp{T) > re- It remains to determine f so that {r{t),v{t)) is part of 
a controlled admissible trajectory of System One first integrates over [0,T] the differential 
equation 

and then take T(t) = (yv{t) + rriit) for t G [0, T]. The final bound f is simply the 

maximum of ||T(t)|| over [0,T]. 

□ 


As a combination of Lemmas 5 and 6, we obtain the following result. 

Corollary 2. For every /d > 0, /i > 0, and initial point {xi,mi) G V~ x there exists a limiting 
value Tmax > 0 depending on {xi,mi) sueh that the following properties hold: 

1) if r > Tmax; the system Sgat is eontrollable for the OIP; and 

2) if T < Tmax; the system Sgat ts not eontrollable for the OIP. 

The limiting value Tmax can be computed by combining a shooting method and a bisection method 
as described in Section 4. 
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3.3 Controllability for DOP 

Let us define a system Sgat associated to Ugat as 




sat • 


'r{t) = v{t), 

' + 0r) 

m{t) = II Tit) II, 


(19) 


where all variables are the same as defined in Eq. (1). For every ^ X x and measurable 
control function r taking values in T(Tmax) with Tmax, we define by r(t, T(t), y^) the corresponding 
trajectory of Egat for some positive times. 

Remark 3. For every eontrolled trajeetory (r(t), m(t)) = T{t,T{t),y^) of the system Egat on 

some finite intervals [0,t/] C Xr with (rf^Vf^mf) = r(t/, T(t/), T/i), the trajeetory {r{t),v{t)^rfi{t)) 
r(t, T(tj—t), rj, —vj, mj) o/Egat runs baekward in time along the trajeetoryV(t^T(t)^yfi) on [Ofif], 
i.e., 

{rit),vit),rhit)) = {ritf -t),mitf -t)) on[0,tf]. (20) 

As a consequence, according to Lemma 6 and Corollary 2, we obtain the following result. 

Corollary 3. For eaeh fi > 0, y > 0, and > 0, given a point (vf^Vf) G V~, there exists a 
^max > 0 depending on (vf^Vf) and rui sueh that the following properties hold: 

1) if T > Tmax; the system Egat is eontrollahle for the eorresponding DOP; and 

2) if T < Tmax; the system Egat is not eontrollahle for the eorresponding DOP. 


4 Numerical Examples 

In this section, we consider two numerical examples, one OIP and one DOP, to compute the limiting 
value Tmax In Corollaries 2 and 3, respectively. The gravitational constant y in system Egat (and/or 
^sat) is 3986000.47 Km^/s^, the radius of the Earth is Te = 6, 374, 000 m, and we consider the 
vertical depth of atmosphere around the Earth is 90,000 m, which means Tc = re + 90, 000 m. 

4.1 A numerical example for OIP 

In order to be able to compute the limiting value Tmax in Corollary 2, we first define the following 
optimal control problem. 

Definition 12 {Optimal eontrol problem (OOP) for OIP ). Given every initial point {xi^mi) G 
V~ X and r > 0, the optimal eontrol problem for OIP eonsists of to steering a satellite by 
t(-) G T(t) on a time interval [0,tj] C Xp along System Egat sueh that, along the eontrolled 
trajeetory {r{t),v{t),m{t)) = V{t,T{t),Xi,mi), the time tf is the first oeeurenee for || r{tf) || = 
rp{r{tf),v{tf)), i.e., || r{t) \\> rp{r{t),v{t)) on [0,t/), and that rp{r{tf),v{tf)) is maximized, i.e., 
the eost funetional is 


Let > 0 be the optimal final time of the OCP for OIP, and let {x{t),m{t)) = T{t,T{t),Xi,mi) 
on [0,t/] be the optimal controlled trajectory with the associated optimal control T{t) G T(t) 
on [0,t/]. One can check, by using Pontryagin Maximum Principle as was done in Ref. [7], that 
II T{t) 11= r on the whole interval [0,tj]. Thus, fixing the initial point {xi,mi) G V~ x we 
have that the final time tf, the trajectory {x{t),m{t)) at each time t G [0, t/], and the final perigee 
distance rp{x{if)) are functions of r. Thus, let us define a function 

s : M+ ^ M, s{t) = rp{x{if)) — re- (22) 
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If one can find Tmax > 0 such that ^(rmax) = 0, then Tmax is the limiting value in Corollary 2. For 
every r > 0, using a shooting method as an inner loop to solve the OCP for OIP, we can obtain 
a value for s(r). Then, using a bisection method as an outer loop, one can obtain Tmax > 0 such 
that 5(Tmax) = 0. Accordiug to Eq.(14) and the objective of the OCP for OIP, place a satellite 
with the initial mass > 0 on a point {vi^Vi) G V~. The optimal controlled trajectory lies 
on a 2-dimensional plane spanned by and Vi. Hence, the limiting value Tmax in Corollary 2 is 
determined only by || Vi ||, || Vi ||, and by the flight path angle rji G [—7r/2,7r/2], i.e., the angle 
between the velocity vector Vi and local horizontal plane, defined by 


r]i = sin 



Assume that a rocket carries a satellite, whose initial mass is rrii = 150 kg^ from the surface of 
the Earth to a point Xi = (vi^Vi) in the unstable region V~ such that || Vi ||= re + 110,000 m, 
II Vi 11= 7879.5 m/s, and rji = 5°. The rocket and the satellite are separated at this point Xi. 
Then, the satellite has to use its own engine to steer itself from the point Xi into the stable region 
. We can see from Eigure 4 that the periodic orbit 73,. has collisions with the surface of the 
atmosphere around the Earth. 



X(m) ,,10^ 


Eigure 4: The periodic orbital 73.. and the optimal controlled trajectory of the OCP for the OIP 
with r = Tmax Starting from Xi. 

We choose the specific impulse of the engine fixed on the satellite as Ap = 2000 s, which implies 
that P = = 5.102 X 10 ^ m ^ where go = 9.8 m^/s. The computed result of the limiting 

value is Tmax = 8.052 N. Thus, in order to be able to insert the satellite from the point Xi into 
stable region the maximum thrust of the engine has to be larger than 8.052 N. The optimal 
controlled trajectory of the corresponding OCP for OIP with r = 8.052 is plotted in Eigure 4 as 
well. To see the numerical results for different initial points, another two points Xi = and 

^2 = are chosen on the periodic orbit 73.. such that 

\\ri 11= re + 379,494 m, || ||= 7,562 m/5, 71 = 4.3517°, 

II r 2 11= re + 599,351 m, || i;2 ||= 7,312 m/5, 72 = 3.0132°. 

Then, the limiting value of Tmax corresponding to the two initial points Xj^ j = 1, 2, are computed 
as 9.037 N and 10.719 N, respectively. We see that the limiting values Tmax are different for 
different initial points on the same periodic orbit. The time history of radius || r{t) || and perigee 
distance rp{x{t)) along the optimal controlled trajectories starting from the initial points Xi and 
Xj, j = 1,2, are plotted in Eigure 5. Since the three points Xi and Xj^ j = 1,2 lie on the same 
periodic orbit 73.., Vp at initial time is the same, as shown in Eigure 5. 
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Figure 5: The profile of || r || and Vp with respect to time for 3 different initial points, i.e., Xi and 
Xj {j = 1,2), for OIP; The subscripts, 7, 1, and 2, correspond to the initial points cci, and X 2 ^ 
respectively. 

4.2 A numerical example for DOP 

A DOP is a powered flight phase of a satellite in the region A x during which a decelerating 
manoeuvre is performed so that the satellite will move to the desired final point Xf = (vf^Vf) G V~ 
at the entry interface (El). The condition at El permits the satellite to have a subsequent safe entry 
flight in atmosphere to a landing site. A typical condition at El, see Ref. [1], is given as: 

II rj 11= tei, II Vf 11= Vei, and rj -vj = VEirEisin{r]Ei), (23) 

where vei = 122,000 m, Vei = 7879.5 m/s, and tjei = —15° denote the norm of position 

vector, the norm of velocity vector, and the flight path angle at El, respectively. In order to 
compute the limiting value Tmax in Corollary 3 for the DOP to a point (vf^Vf) in V~ ^ we first 
define the below optimal control problem. 

Definition 13 {Optimal control problem (OCP) for DOP). Given every final point x f = {vf^Vf) G 
V~ and r > 0; let mi > 0 be the initial mass of a satellite, the optimal control problem for DOP 
consists of steering the satellite by t(-) G T(r) on a time interval [0,tf] C Xr subject to the system 
Esat such that, along the controlled trajectory {r{t),v{t),m{t)) = T{t,T{t),r f, —Vf,mf) (mf > 0 is 
free) of System Egat; the time tf is the first occurence for || r{tf) \\ = rp{r{tf),v{tf)), mi = m{tf), 
and rp{r{tf),v{tf)) is maximized, i.e., the cost functional is the same as Eq.(21). 

Given every initial mass m^ > 0 and final point {rf,Vf) in V~, let 7/ > 0 be the optimal final 
time of the OCP for DOP, and let {x{t),fh{t)) = V{t,T{t),r f, —Vf,mf) be the optimal controlled 
trajectory associated to the control t{t) G T(t) on [0,7/]. Then, the same as the OCP for OIP, 
the perigee distance rp{x{tf)) is a function of r. Let us define a function 

s : M+ ^ M, s{t) = rp{x{if)) - rc. (24) 

Then, according to Eq. (20), in order to compute the limiting value r^ax in Corollary 3, it suffices 
to combine a shooting method and a bisection method to compute the value Tmax such that 

'^('^max) — 0. 

What we developed in this paper is applicable not only for low-thrust control systems but also 
for high-thrust control systems if only the thrust is finite instead of impulsive. Thus, we consider the 
space shuttle’s parameters in Refs. [3, 10]. The initial mass is 95,254.38 kg. The specific impulse 
of the engine is 313 s that means p = 3.26 x 10“^. The numerical result is Tmax = 14,004.62 N. 
Note that the propulsion for a space shuttle is provided by the orbital manoeuvring system (CMS) 
engines, which produce a total vacuum thrust of 53,378.6 N, see Refs. [3, 10]. Thus, according to 
Lemma 5, for every initial point Xi in the space shuttle can reach the El condition in Eq.(23) 
by admissible controlled trajectories of the system Egat if the satellite takes enough fuel. The 
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periodic trajectory and associated optimal controlled trajectory with r = r^ax are illustrated 
in Figure 6. The profile of || r || and Vp along optimal controlled trajectories for the DOP with 


X 10® 



Figure 6: The periodic trajectory determined by the El condition in Eq.(23) and the optimal 
controlled trajectory of the OCP for DOP with r = Tmax- 

T = Tmax, Tmax + 100 N^ and Tmax “ 100 N^ are illustrated in Eigure 7. We can see from Eigure 7 
that the optimal controlled trajectory of the OCP for DOP with t = Tmax +100 N is an admissible 
controlled trajectory in A and the final point lies in . While, the optimal controlled trajectory 
of the OCP for DOP with t = Tmax — 100 N cannot reach a point in by admissible controlled 
trajectories of the system Egat- 



Eigure 7: The profile of r =|| r || and Vp along the optimal controlled trajectory of the OCP for 
the DOP with t = Tmax, Tmax + 100 N, and Tmax “ 100 N. 


5 Conclusion 

The controllability property of the Keplerian motion around the Earth in the periodic region V is 
established in this paper. According to the state constraint that the radius of the Keplerian motion 
has to be larger than the radius of the surface of atmosphere around the Earth, the periodic region 
is separated into two sets: and V~. The controlled motion in the set is the typical OTP 

and we obtain that the motion is controllable in the set for any positive maximum thrust. 
Moreover, we obtain that there exists a limiting value of Tmax > 0 depending on initial point (final 
point, respectively) such that the Keplerian motions for OIP (DOP, respectively) is controllable if 
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r > Tmax- Finally, two numerical examples are simulated to show that a shooting method and a 
bisection method can be combined to compute the limiting value for the bound on the thrust. 


6 Appendix 

In this section, we provide two sets of coordinates for points in the periodic region V (see Ref. 
[8, 15] for all the results given here). 

Definition 14 {Classical orbital elements (COE)). For x G V, define the following functions: 


a{x) 

II 

1 

(25) 

i{x) 

- cos- IIII 

II ^ mil. II ’ 

(26) 

uj{x) 

- cos- II ^"-"11 

Hi nil nil’ 

(27) 

Tt{x) 

- cos- ll'F-ll 

111. III! nil’ 

(28) 


where n = Iz x h with Iz = [0,0,1]^. The quantity a{x) is called the semi¬ 

major axis of the orbit whose shape is thus determined by a{x) and e{x). The angles i{x), 
uj{x) and ft{x) are called the inclination of the orbit the argument of perigee of the orbit jx 
and the right ascension of the ascending node of the orbit ^x respectively. Then, the variables 
{a{x),e{x)fi{x),uj{x),Vt{x),0{x)) are called the classical orbital elements of the orbit ^x- 

(Note that the set of COEs is singular if e = 0 and i = 0 , tt.) 

Definition 15 {Modified equinoctial orbital elements (MEOE)). For x G V, define the following 
functions: 


P{x) 

= a{x){l- e{xf)/fi, 

(29) 

(^) 

= e{x) cos{uj{x) P fl{x)), 

(30) 

ey{x) 

= e{x) sm{uj{x) P Q{x)), 

(31) 

hx (^) 

= tan(i(cc)/2) cos(r^(cc)). 

(32) 

hy{x) 

= tan(i(cc)/2) sin(f2(x)). 

(33) 

fix) 

= uj{x) P ft{x) P 0{x), 

(34) 


where {a{x),e{x)fi{x),uj{x),Vt{x),0{x)) are the COE defined previously. Then the 6-tuple z = 
{P^Cx^Cy^hx^hyfi) G X S gathers the so-called modified equinoctial orbit elements (MEOE), 
Moreover, we also have that 


r 


V 


P 


(1 + — hy) cos I + 2hxhy sin I 

(1 - + hy) sin / + 2hxhy cos / 

2hx sin / — 2hy cos/ 


C 


2hxhy{ex + cos/) - (1 + - hy){ey + sin/) 

-2hxhy{ey + sin/) + (1 - + h‘l){ex + cos/) , 

2hx{ex + cos/) + 2hy{ey + sin/) 


(35) 

(36) 


where C = 1 + and W = 1 + cos / + sin/. Note that e = y and P = / pi. 

Thus, let us define the set 

Z = e X S : P > 0 and 0 < < l}, 

then the transformation {r,v) : Z ^ V, z {r{z),v{z)) is a covering map. Hence V is arc- 
connected if Z is. is sufficient 
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